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We investigate the generation of gravitational waves in the hybrid quintessential inflationary 
model. The full gravitational-wave energy spectrum is calculated using the method of continuous 
Bogoliubov coefficients. The post-inflationary kination period, characteristic of quintessential infla- 
tionary models, leaves a clear signature on the spectrum, namely, a peak at high frequencies. The 
maximum of the peak is firmly located at the MHz-GHz region of the spectrum and corresponds to 
f^Gw — 10~^^. This peak is substantially smaller than the one appearing in the gravitational-wave 
energy spectrum of the original quintessential inflationary model, therefore avoiding any conflict 
with the nucleosynthesis constraint on ficw. 
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I. INTRODUCTION 

Gravitational waves of cosmological origin are at 
present the object of an important research effort, as 
people realize that they will provide us with a unique 
telescope to the very early Universe, giving us informa- 
tion not easily available by any other means. Indeed, it 
would be of the utmost importance to be able to find di- 
rect signals coming from these earliest of times, in order 
to improve our understanding about inflation, preheat- 
ing and reheating mechanisms, post-inflationary phase 
transitions, topological defects of grand unified theories, 
and string theory, among other issues relevant for both 
cosmology and high-energy physics. 

In this work we investigate gravitational-wave genera- 
tion within a model recently proposed by Bastero-Gil et 
al. [H , a modification of the original quintessential infla- 
tionary model of Peebles and Vilenkin 0] . In the original 
model, a unified description of infiation and dark energy 
is achieved with a single scalar field playing both roles 
of infiaton and quintessence. In such quintessential infla- 
tionary models reheating does not proceed in the usual 
way, through the complete decay of the scalar field into 
quanta of other fields. In fact, since the potential of the 
scalar field has no minimum, infiation is followed not by 
coherent oscillations of the scalar field, but rather by a 
kination period Q , during which the evolution is dom- 
inated by the kinetic energy of the scalar field ^, which 
behaves approximately as stiff matter. Reheating pro- 
ceeds then by gravitational particle production taking 
place at the transition from the inflationary to the kina- 
tion period Q. This mechanism is, however, quite ineffi- 
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cient and may lead to cosmological problems associated 
with large isocurvature fluctuations and overproduction 
of gravitinos and moduli fields In order to recover 
the usual reheating mechanism, Bastero-Gil et al. [l| in- 
troduce another scalar field Xi coupled to the original in- 
fiaton/quintesscncc field with a hybrid- like potential 
0. In such a hybrid quintessential inflationary model, 
during the kination period driven by the kinetic energy 
of the field 0, the field x oscillates around the minimum 
of the hybrid-like potential, completely decaying into rel- 
ativistic particles, thus reheating the Universe. Since this 
reheating mechanism is more efficient than gravitational 
particle production, the radiation-dominated epoch of ex- 
pansion of the Universe begins earlier and, therefore, the 
kination period is shorter in this hybrid model than in 
the original one. 

It is known that a phase of evolution in the early Uni- 
verse with an equation of state stiffcr than radiation leads 
to a sharp increase of the gravitational- wave spectral en- 
ergy density parameter Q.^^ in the high-frequency region 
of the spectrum [sl-fioj. This peak in the gravitational- 
wave spectrum is unavoidable in quintessential infiation- 
ary models, in which inflation is followed by a period 
of kination In the original quintessential infla- 

tionary model of Peebles and Vilenkin the height 
of the peak depends on the number of scalar degrees 
of freedom, N^, whose decay into fermions triggers the 
onset of a gravitational reheating of the Universe. For 
the minimum allowed value of Ng, of the order of 10^, 
the gravitational-wave spectral en ergy density parame- 
ter l](3w could be as high as 10~^ [ll|, several orders of 
magnitude higher than in standard infiation. Clearly, in 
the hybrid quintessential inflationary model, in which the 
kination period is shorter than in the original model, one 
expects the high-frequency peak of the spectrum to be 
smaller. 

Our main aim in this paper will be to calculate, within 
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the hybrid quintessential inflationary model, the present- 
day spectrum of the cosmological gravitational waves 
generated during the evolution of the Universe and to 
show that its careful analysis will allow us to extract 
much information about the properties of the potential 
driving both the primordial inflationary stage and the 
present-day accelerated expansion. To take into account 
the expansion of the Universe during the transitions be- 
tween the different stages of its evolution, characterized 
by different equations of state, we use throughout the 
method first developed by Parker, in his seminal paper 
of 1969 [iHl- This method, applied by Parker to particle 
production in an expanding Universe and then extended 
to the case of gravitons jl6l - ll8l |. is based on the time 
evolution of the Bogoliubov coefficients, which obey ap- 
propriate differential equations. The numerical integra- 
tion of these equations will immediately allow us to con- 
struct the full spectrum of the gravitational waves (see 
Refs. [l9l - l22j for applications of this method to several 
cosmological models). 

Our paper is organized as follows. In the next section 
we describe the hybrid quintessential inflationary model 
of Bastcro-Gil et al. [l[ and write the equations of motion 
for the different stages of the expansion. We also describe 
the simple phenomenological mechanism responsible for 
the reheating of the Universe, acting during the kination 
period that follows inflation. In section IIII[ we present 
the aforementioned method of the continuous Bogoliubov 
coefficients to calculate the spectrum of the gravitational 
waves, integrate the corresponding equations and com- 
pute the gravitational-wave spectral energy density pa- 
rameter for frequencies ranging from about 10""'^^ rad/s 
to about 10^° rad/s. The influence on the spectrum of 
the various parameters defining the hybrid quintessen- 
tial inflationary potential and the resulting constraints 
are then carefully analyzed in the different frequency re- 
gions. In section HVl we summarize our main conclusions. 

II. HYBRID QUINTESSENTIAL INFLATION 

A. The hybrid potential 

As shown by Peebles and Vilcnkin 0], the unifica- 
tion of inflation and dark energy within a single frame- 
work, which they termed quintessential inflation, can be 
achieved with a single scalar fleld with potential given by 

r A^(</)4 + m4), for (/XO, 

= <^ (1) 
[ A0M8(04 + m4)-i, for (/) > 0, 

where and M are constants. The value of A^ is con- 
strained by recent measurements of the cosmic microwave 
background and large-scale structure to be of the order of 
10"^'^ (see Ref. [2^ for a derivation of such constraints for 
different models of inflation) . On the other hand, agree- 
ment with the measured value of today's dark-energy 
density [g^] requires M to be of the order of lO^^^rrip. 
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FIG. 1: Potential U{4>,x) for m = 0.004 rrip, = 10"", 
M = 1.7 X 10~" mp. Ax = 1 and g = 0.025. 

Recently, Bastero-Gil et al. [l|] proposed a modiflca- 
tion of the original quintessential inflationary model, in 
which the reheating of the Universe is driven not by gra- 
vitational particle production but rather by the simpler 
mechanism of the decay of an oscillating massive fleld. 
This is achieved by the introduction of an extra scalar 
fleld X with a hybrid-like potential given by 

C/(0, x) - V{4>) + \gW4>' - m^) + \x^x\ (2) 

where V{4>) is given by Eq. ([1]) and g, m, and A^ are 
some parameters to be specifled. Note that for |^| ^ 
m the potential ^ has just one minimum at x = 0, 
while for \(f)\ < m it has two minima located aX x = 
±^g'^{m? -4)'^)l\ (see Fig.p. 

In our analysis of this hybrid quintessential inflationary 
model we will divide the evolution of the Universe into 
four stages. 

In the first stage of evolution, ^ — m and x is as- 
sumed to be at rest at the minimum of the potential 
(x = 0)- Initial conditions are such that the potential 
energy of the scalar field cj) dominates the evolution of 
the Universe, yielding a period of chaotic inflation. By 
the end of this stage of evolution, the kinetic term in the 
energy density of (p becomes dominant, giving rise to a 
period of kination, during which the scalar field behaves 
as stiff matter with equation of state p = wp, w = 1. 

During the second stage of evolution, which takes place 
for 101 < TO, the origin of the potential becomes unsta- 
ble for the scalar field x- As a result, this field rolls 
down toward one of the new minima of the potential, lo- 
cated at X = -^\/ g'^im? — (/)^)/A^, and starts to oscillate 
around it. Of course, this behavior of the scalar field 
X is expected only if the duration of the second stage 
of evolution is much greater than the characteristic re- 
sponse time of the field to changes in the potential. As 
shown in Ref. [l[, this requirement can be formulated as 
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a constraint on the values of the parameters g and to, 
namely, 



(3) 



We also demand that the motion of the scalar field 
is not influenced significantly by i-S-, we impose the 
condition that the energy density of x is much smaller 
than the kinetic term in the energy density of (j) during 
the second stage of evolution. This condition can be also 
translated into a constraint on the values of the param- 
eters appearing in the potential ^ , namely [l| , 
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Note that this last condition implies that kination ex- 
tends throughout the second stage of evolution. 

The third stage of evolution takes place for ^ to. 
For such values of the scalar field (j), the potential ([2]) 
has again only one minimum in the x direction, located 
at the origin. The scalar field x oscillates around this 
stable minimum, transferring its kinetic energy, acquired 
during the second stage of evolution, to a radiation fluid, 
thus reheating the Universe. This decay of the oscillating 
scalar field x into radiation is achieved by the introduc- 
tion of a phenomenological dissipative coupling, between 
the scalar field y and the radiation fluid, proportional to 
the mass of x H Ull , 



(5) 



fi being the proportionality constant. The above- 
mentioned condition that the motion of the scalar field 
(j) is not influenced signiflcantly by ^ also imposes a con- 
straint on the values of fj., namely 



X a 

i. 
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Note that this constraint, as well as constraints ([3]) and 
(|4|), was derived assuming that to <IC nip^. During the 
third stage of evolution, the energy density of x de- 
cays away rapidly. Gradually, as the energy density of 
radiation increases, kination gives place to a radiation- 
dominated Universe. 

The fourth stage of evolution extends from the begin- 
ning of the radiation-dominated era to the present epoch. 



^ The constraints l[3]l, l|4]l, and (|6]l differ from the ones derived in 
Ref. Q by some numerical factors. This is due to the fact that 
Bastero-Gil et al. estimated the scale-factor growth between the 
end of the inflationary period and the beginning of the second 
stage of evolution to be of the order of 8, while our numerical 
simulations show this growth to be half of this value. We take 
into account this factor 2, since the constraints involve third and 
sixth powers of the scale-factor growth. 



The scalar field x is assumed to have decayed away com- 
pletely in the previous stage of evolution. A matter com- 
ponent (dark and usual, baryonic, matter) is introduced 
into the equations of motion, giving rise to a interme- 
diate matter-dominated period in the evolution of the 
Universe. The scalar field (p, which in the first stage of 
evolution played the role of inflaton, stays now practically 
constant and, at late times, begins to dominate the evo- 
lution of the Universe, giving rise to the present epoch of 
accelerated expansion. Therefore, the scalar field cf> plays, 
within this model, both roles of infiaton and quintessence. 



B. The evolution of the Universe 

Let us now present the equations of motion for the 
different stages of evolution. 

For the first and second stages of evolution, the equa- 
tions of motion are 
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where we have assumed a flat Friedmann-Robertson- 
Walker metric, a is the scale factor, TOp is the Planck 
mass^, and a dot denotes a derivative with respect to the 
cosmic time t. 

A Runge-Kutta method is used to solve the system 
of differential equations ([Zl)-®, while Eq. (fTO|) is used 
as a constraint equation to check the accuracy of the 
numerical solution. 

As initial condition for the scalar field <f> we choose 
~ —5 TOp, which guarantees enough inflation. The 
scalar field x is located near the origin; for numerical 
convenience, we assume Xii = Xii = 0- The initial val- 
ues of t he other vari ables are chosen to be an = 1 and 
4>ii = y^2U{4>ii,Xii), while dn is fixed by the Friedmann 
equation ([T0| . As mentioned above, the parameters 
and M of the potential V{(j)) are constrained by observa- 
tions to be of the order of 10~^^ and 10~^^TOp, respec- 
tively. The first stage of evolution ends when </> = — to; 
for the parameter to we choose to < 10^^ TOp. 

For such values of the initial conditions and of the 
parameters A^, M, and to, the first stage of evolution 



^ In the first, second, and third stages of evolution we use the na- 
tural system of units, with h = c = 1 and rrtp = G"^/^ = 1.22 X 
10^^ GeV, while in the fourth stage we use the International 
Systems of Units. 
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FIG. 2: Time evolution of the equation-of-state parameter 
w = p/ p, from the beginning of the inflationary period until 
the present epoch, for = 10~^^, M = 1.7 x 10~^*mp, 
g = 0.025, m = 0.004 mp, = 1, and p. = 0.1. In this 
figure, instead of the cosmic time t, we use the variable u = 
— In(l-l-z), where z is the redshift. The first stage of evolution 
corresponds to —139 < u < —63, the second stage to u « —63, 
the third stage to —63 u < —60, and the fourth stage to 
—60 < It ^ 0. The epochs of inflation, kination, radiation 
domination, matter domination, and dark-energy domination 
can be clearly identified in this figure. 
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FIG. 3: Allowed values of the parameters m and g (shaded 
region), for A^ — 10"^'^ and A;^ — 1. 



lasts for about \(f and the scale factor grows about 
32 orders of magnitude. The equation-of-state parame- 
ter w = pI p changes gradually from w = — 1 (inflation) 
to w = 1 (kination), as shown in Fig. [5] 

At the beginning of the second stage of evolution, the 
origin becomes unstable for the scalar field x- We per- 
turb slightly this field, making it roll towards the tempo- 
rary minimum located at Xmin(0) = +\/g^(?Ti^ — 0^)/A-^. 
More exactly, we choose Xi2 = O.OlXmin(O) and Xi2 = 0. 
The parameters g, to, and A^^- should satisfy the con- 
straints ^ and (|H). We also demand that 5 ^ 1, < 1, 
and TO < lO^^JTip. If, for instance, we choose \^ ~ 1, 
then the values of g and to are restricted to the shaded 
region shown in Fig. |31 For such values of the parame- 
ters the second stage of evolution lasts typically for about 
lO^tp (see Fig. 15). 



FIG. 4: Time evolution of the scalar field x during the sec- 
ond and third stages of evolution, for the same values of the 
parameters as in Fig. [2] The second stage of evolution lasts 
about 7x 10^ tp; the third stage of evolution, which lasts about 
5 X 10^° fp, is only partially shown is this figure. 



During the third stage, the evolution of the Universe 
is described by the set of differential equations 
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(11) 

(12) 
(13) 
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(15) 



where pr is the energy density of radiation and F^ is the 
dissipative coefficient given by Eq. ([5|). 

Again, we use a Runge-Kutta method to solve the sys- 
tem of differential equations (|lip - p^ . Equation ([15]) 
is used to check the accuracy of the numerical solution. 
Since any pre-existing radiation fluid would have been 
diluted during inflation, we choose the energy density of 
radiation at the beginning of the third stage of evolution 
to be zero, pr,i3 = 0. The dissipation parameter p should 
satisfy the constraint ([5]). For A^ = 1, to = 0.004 mp, 
and g = 0.025 this implies p ^ 10~^. 

As already mentioned above, during the third stage of 
evolution, the scalar field x oscillates around the mini- 
mum of the potential located at x = (see Fig.|3]), trans- 
ferring its energy into the radiation fluid (see Fig.[S]), thus 
reheating the Universe. The reheating temperature de- 
pends on the values of the parameters of the model [l[ . 
For the above-mentioned values of the parameters, our 
numerical simulations show that the reheating tempera- 
ture is of the order of 10^* GeV. 

For low values of the parameter /i, the scalar field x os- 
cillates for a long period of time before its energy is com- 
pletely transferred to the radiation fluid. In such cases. 
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the previous stage of evolution. Therefore, Eq. (|14l) . 
or Eq. (PH)) . can be integrated exactly, yielding ~ 

Pr,o(ao/a)^j where pr,o — 4.13 X 10'" J/ni and flo are, 
respectively, today's values of the energy density of radi- 
ation and of the scale factor. In this stage of evolution a 
pressureless matter component is introduced, which ac- 
counts for the usual (baryonic) matter and also for dark 
matter, with energy density pm = Pm,o(ao/a)'^j where 
Pmfi — 2.34 X 10~^° 3/tc? is today's value of the energy 
density of matter. Taking all this into account, the equa- 
tions of motion for the fourth stage of evolution are given 
by 



FIG. 5: Time evolution of the energy densities of radiation, <l>, 
and X at the beginning of the third stage of evolution (which 
lasts about 5 x 10^" tp), for the same values of the parame- 
ters as in Fig. [5] Energy is transferred from the scalar field 
X to radiation, while the evolution of the Universe is domi- 
nated by the scalar field (j). Later in this stage, the energy 
density of radiation becomes dominant; kination gives place 
to a radiation-dominated Universe. 



because of the difficulty to accurately follow many oscil- 
lations of X with our numerical code, a different strat- 
egy has to be used. We start by numerically solving 
Eqs. pT|) - ([T^ . After a short period of integration, the 
terms g^m^x? 1'^ £ind A^x"'/4 in the potential J7((/), x) be- 
come much smaller than ^■^ 0^x^/2, implying that the en- 
ergy density of the scalar field x is given, within a good 
approximation, by 



~^ I 2 i2 2 



(16) 



Taking into account that remains practically constant 
during an oscillation of x^ the oscillations of the latter 
field can be considered to be of the simple harmonic type, 
for which p^ = (x^) — 5^0^(x^)j where (...) denotes the 
average over one oscillation. Using this approximation, 
Eqs. ([Tl]) -(fT5 |) can be re-written as 
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,(22) 
(23) 

(24) 



where G is the gravitational constant and c is the speed 
of light. 

This set of differential equations describes a stage of 
evolution in which the Universe is dominated, consecu- 
tively, by radiation, matter and the inflaton/quintessence 
field (j). As already mentioned above, in order to obtain 
agreement with observations [l^j, the parameter M of 
the potential V{4>) has to be chosen such that 



(25) 



where pde,o — 6.20 x 10~ J/m is today's value of the 
energy density of dark energy. Within the present model, 
the value of AI is typically of the order of 10~^^ rrip. 

Note that the values of pde.O! Pm,o, and pr.o used in 
this article imply a value of the Hubble constant of Hq = 
71 kms"iMpc"^ 



V+^PX+ Pr 
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Pr + 4-pr = r^Px' 
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We then solve these equations, instead of Eqs. (fTT|) -([T4 l) . 
until the end of the third stage of evolution. 

Let us now turn to the fourth stage of evolution. The 
scalar field x has already decayed away completely in 



III. GRAVITATIONAL WAVES 

Gravitational waves are generated in an expanding 
Universe, giving rise to a spectrum extending over a wide 
range of frequencies, from about 10""'^^ rad/s to about 
10^° rad/s (26H32. 

In this article, we calculate this spectrum for the hy- 
brid quintessential inflationary model using the method 
of the continuous Bogoliubov coefficients. This method 
can be summarized as follows (for details, see Refs. [l6l - 
[2^). The number of gravitons at a certain moment of the 
evolution of the Universe is given by the squared Bogoli- 
ubov coefficient, f^i^ = {X — 1^)^/4, where the functions 
X{t) and Y{t) are solutions of the system of differential 
equations 



X = -toJo-Y, 
a 



(26) 
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ao and luq are today's values of the scale factor and the 
gravitational-wave angular frequency, respectively. The 
above system of equations is integrated with initial con- 
ditions X{tii) = Y{tii) = 1, corresponding to the ab- 
sence of gravitons at the beginning of the first stage of 
evolution. The scale factor a{t) and its first and sec- 
ond derivatives d{t) and d{t) are determined from the 
evolutionary equations presented in the previous section. 
Knowing /3'^{t), we can compute the gravitational- wave 
spectral energy density parameter, flaw, which is defined 
as 



,4o2 



(28) 



where H{t) is the Hubble parameter. Evaluating all the 
quantities in the above expression at the present time, for 
angular frequencies ranging from 1.4 x 10~^^ rad/s (cor- 
responding to a wavelength equal, today, to the Hubble 
distance) to about 10^° rad/s (corresponding to a wave- 
length equal to the Hubble distance at the end of the in- 
flationary period), yields the gravitational- wave energy 
spectrum. It is worth emphasizing that this spectrum 
is fully determined by the evolution of the scale factor 
from the beginning of the inflationary period until the 
present time. A change in the behavior of the scale fac- 
tor, due to a change of the equation of state of the Uni- 
verse (see Fig. ^ , leads to a modification of the slope of 
the gravitational-wave spectrum at a certain frequency. 

Let us now compute the gravitational-wave spectrum 
of the hybrid quintessential inflationary model using 
the formalism of the continuous Bogoliubov coefficients. 
First, we have to specify the values of and M, satisfy- 
ing the observational constraints on inflation and dark 
energy, and the values of g, m, Aj^-. and fi, satisfy- 
ing the constraints ([S])-®. In Fig. IH] we show the full 
gravitational-wave spectrum for the hybrid quintessen- 
tial inflationary model for two different values of the pa- 
rameter g, corresponding to the minimum and maximum 
values allowed by the constraints (jH]) and ^ for fixed 
values of A^, m, and A^. 

As we can see from Fig. [51 the spectrum is naturally 
divided in three regions. 

In the low- frequency region, luq ranges from 1.4 x 
10^^'' rad/s to about 2TTH{t„i)a{tjn)/ao ~ 10"^^ rad/s, 
where the upper limit corresponds to today's value of 
the angular frequency of a gravitational wave which had 
a wavelength equal to the Hubble distance at the time 
trn when the energy density of radiation became equal 
to the energy density of matter. In this region of the 
spectrum, the gravitational- wave spectral energy density 
parameter ilcw rises rapidly as ujq decreases, due to an 
extra graviton production during the transition between 
the radiation- and the matter-dominated eras and the 
subsequent matter- and dark energy-dominated eras. For 
Wo = 1.4x 10~^^ rad/s the spectrum satisfies a constraint 
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Gravitational-wave spectrum for 
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to the minimum and maximum values of g allowed by the 
constraints Q and Q, namely, g — 0.004 and g = 0.025. In 
both cases AI is of the order of 10^^* m„. 



derived from measurements of the cosmic microwave 
background radiation, namely, ficw < 1-4 x 10~^^ [s^ . 

In the intermediate-frequency region, the angu- 
lar frequency ranges from about 10^^^ rad/s to 
2TTH{tr)a(tr)/ao, where tr corresponds to the transi- 
tion between the kination and radiation-dominated eras, 
occurring at the end of the third stage of evolution. 
Depending on the duration of the kination era, the 
intermediate-frequency region extends up to values of 
about (lO'^ — 10^) rad/s. Note that the spectrum has, in 
this region, a nonconstant slope, due to the fact that in- 
flation is quasi-exponential'^ 0, HI] . Several bounds on 
the gravitational-wave spectral energy density parame- 
ter should be satisfied in this region of the spectrum, 
namely, from timing observations of millisecond pulsars, 
flew < 4 X 10"^ for Wo = 2.5 x 10"^ rad/s from 
Doppler tracking of the Cassini spacecraft, flew < 0.028 
for Wo = 7.5 X 10~^ rad/s [IHl, and from the Laser In- 
terferometer Gravitational Wave Observator y ( LIGO), 
r^Gw < 6.9 X 10-6 for Wo « (10^ - 10^) rad/s j3|. 

In the high-frequency region of the spectrum, ranging 
from (10"^ — 10^) rad/s to about 10^° rad/s. one observes a 
sharp rise of the gravitational-wave spectral energy den- 
sity parameter Qqw- This rather interesting feature of 
the spectrum is due to the existence of a kination pe- 
riod (with equation of state p = wp, w = 1) between 
the end of inflation and the beginning of the radiation- 
dominated era, whose duration determines the height of 
the peak (Toj . 

Our numerical simulations show that the duration of 
the kination period increases as the value of the param- 
eter g decreases (see Fig. [7|) . This can be easily under- 



^ In exponential inflation (cosmological constant) the slope of the 
spectrum in the intermediate-frequency region is zero [27ll, while 
in power- law inflation the slope is constant and negative [3011 . 
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FIG. 7: Time evolution of the equation-of-state parameter 
w = p/p for g = 0.004 and g = 0.025. The duration of the 
kination period in the former case is much longer, implying 
a higher peak in the gravitational-wave spectrum. In this 
figure the values of the parameters are the same as in Fig. |6] 
and u = — ln(l -|- z), where z is the redshift. 



stood as follows. If one decreases g (for fixed values of m 
and X^), the value of x at the minimum of the potential 

also decreases (recall that Xmin = ± \J Q^ijn^ — 'P)/^x ) ' 
meaning that less energy is acquired by this field as it 
oscillates around the minimum during the second stage 
of evolution. As a consequence, less energy is available 
to be transferred, during the third stage of evolution, 
from the scalar field x to the radiation fluid, leading to a 
lower reheating temperature. This, in turn, implies that 
after the complete decay of x a longer period of time is 
required for the energy density of radiation to become 
greater than the kinetic energy of the scalar field 0, i.e., 
the beginning of the radiation-dominated era is delayed 
and. consequently, the kination period becomes longer. 
In short, we could say that reheating becomes less effi- 
cient as g decreases. 

As we have seen in Sect. HH g is not a free param- 
eter, its value is bounded from above by the condition 
that the motion of cf) is not affected by x [see Eq. (jl|)] 
and from below by the condition that the scalar field 
X responds quickly enough to changes in the potential 
t/(0, x), rolling down toward one of the new minima and 
oscillating around it [see Eq. ([3])]. Therefore, within the 
hybrid quintessential inflationary model, the duration of 
the kination period and, consequently, the height of the 
peak in the gravitational- wave spectrum, cannot be freely 
adjusted, they are limited by the allowed values of the 
parameter g. 

What in the previous two paragraphs was said rela- 
tively to the dependence of the gravitational-wave spec- 
trum on the value of g could also be said, with the neces- 
sary adaptations, about the parameters m and A^. The 
dependence of the spectrum on the parameter m is il- 
lustrated in Fig. m There, two spectra are shown, cor- 
responding to the minimum and maximum values of m 
allowed by the constraints ^ and Q for flxed values 
of A0, g, and Aj^,. For both values of m, the parameter 



FIG. 8: Gravitational- wave spectrum for = 10^ , 5 = 
0.01, A;^ = 1, and ^ = 1. The two curves correspond 
to the minimum and maximum values of m allowed by the 
constraints ^ and (U), namely, m = 2.5 x 10~'^ mp and 
m — 0.01 TTip. In both cases M is of the order of 10^^* m^. 
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FIG. 9: Gravitational-wave spectrum for A^ = 10~ , g = 
6.3 X 10^*, m = 0.01 mp, and ^ — 0.1. The three curves 
correspond to Ax = 10^^, 10~^, 1. The minimum value of 
A;^, allowed by constraints ((4)) and ([6]), is 1.5 X 10"^ In all 
cases M is of the order of 10"^'* nip. 



/i is chosen such that constraint ^ is satisfled. As we 
see in this flgure, the height of the peak, located at high 
frequencies, increases as the parameter m decreases. In 
Fig. [HI three gravitational- wave spectra are shown for dif- 
ferent values of A^^, and fixed values of A^, g, m and /x. 
As expected, the height of the peak increases as the pa- 
rameter A^ increases. In all cases considered above, the 
parameter M is of the order of 10~^^mp. 

Let us now turn to the analysis of the influence of 
the dissipation parameter /i on the duration of the ki- 
nation period and, consequently, on the height of the 
high-frequency peak of the gravitational-wave spectrum. 
Remember that the parameter ^ is bounded from below 
by the condition that the motion of (fi is not affected by 
X [see Eq. @]. We also require that this parameter is 
smaller than a critical value, ^crit, above which the os- 
cillatory motion of the scalar field x during the third 
stage of evolution would become over-damped. For val- 
ues of the dissipation parameter ^ lying in the interval 
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FIG. 10: Gravitational- wave spectrum for — 10~^^, g — 
6.3 X 10~*, m — 0.01 nip, and Ax = 1. The three curves 
correspond to n = Pcrit = 1.93, fi = 10^^, and /i — 10~®. 
Constraint ([5} requires that S> 10^^. In all cases M is of 
the order of 10"^'' rrip. 



0.2 A^^ {m^/mf^^g^^^X-^ < < /icrit, the duration of 
the kination period decreases as /i decreases. This is due 
to the fact that for smaUer values of n the energy transfer 
of the scalar field x to the radiation fluid proceeds slower, 
with the consequence that the energy density of radiation 
begins to decrease as a~'* later. Therefore, after the com- 
plete decay of the field X: the energy density of radiation 
is higher in the case of small fi and the time it takes 
for this energy density to dominate the dynamics of the 
Universe is shorter, implying that the kination period is 
also shorter. In Fig. [10] we plot three gravitational- wave 
spectra for different values of /x, holding fixed the other 
parameters of the model. We see that the height of the 
peak in the high-frequency region of the spectrum in- 
creases as the dissipation parameter increases, achieving 
its maximum value for fi = /id-it- 

To complete our analysis of the influence of the various 
parameters of the potential U{(j), x) on the gravitational- 
wave spectrum, we point out that the parameter de- 
termines the overall vertical displacement of the spec- 
trum (see Fig. [TT|) . Note that the value of A^ also affects 
the height of the peak, contrarily to the situation in origi- 
nal quintessential inflationary model, where only the low- 
and intermediate-frequency regions of the spectrum were 
pushed down (up ) by an decrease (increase) of A^ (see 
Fig. 4 of Ref. 01). 

To finish this section, let us point out that, for the 
allowed values of the different parameters of the model, 
all gravitational-wave spectra satisfy, by far, the integral 
bound imposed by primordial nucleosynthesis js^ l37l|. 

oo 

J ricwdlnu} < 1.1 X 10-^ (29) 



where a;„ « 10 ^ rad/s. 



-15 -10 -5 5 10 

log [(Jo/ (rad. s"^ ) ] 

FIG. 11: Gravitational-wave spectrum for g = 0.002, m = 
0.01 TTip, Ax = 1, and fi = 0.1. The three curves correspond 
to A^ = 10"^"*, 10"", 10"^^. In all cases M is of the order of 



IV. CONCLUSIONS 

In this work we have investigated the generation of 
gravitational waves in the hybrid quintessential infiation- 
ary model. In this model, recently proposed by Bastero- 
Gil et al. [l|, a unified description of inflation and dark 
energy is achieved with a single scalar field playing both 
roles of inflaton and quintessence, while reheating takes 
place in the usual way, through the complete decay of an- 
other scalar field, Xi rather than by gravitational particle 
production, as in the original quintessential infiationary 
model of Peebles and Vilenkin [2[ . 

The hybrid-like quintessential inflationary potential, 
U((j),x)^ contains several parameters which are con- 
strained by measurements of the cosmic microwave back- 
ground radiation and dark energy, as well as by certain 
restrictions imposed on the behavior of the scalar field Xj 
namely, that it responds quickly enough to changes in the 
potential U (0, x) and that it docs not influence signifi- 
cantly the evolution of the inflaton/quintcsscnce field 0. 

Another relevant parameter of the model is the dissi- 
pation parameter which regulates the rate at which en- 
ergy is transferred from the scalar fleld x to radiation. 
This parameter is constrained from below by the condi- 
tion that the motion of (f> is not influenced signiflcantly 
by X and from above by the requirement that the oscilla- 
tions of X around the minimum of the potential are not 
over-damped. 

For the allowed values of the above-mentioned param- 
eters, we have calculated the full gravitational-wave en- 
ergy spectrum using the method of continuous Bogoli- 
ubov coefficients. Such spectra carry, at high frequen- 
cies, a clear signature of quintessential infiation, namely, 
a sharp rise of the gravitational- wave spectral energy den- 
sity parameter Qcw This distinctive feature of the spec- 
trum is due to the existence of a post-inflationary kina- 
tion period, during which the kinetic energy of the in- 
flaton fleld dominates the evolution of the Universe . 
As we have shown in this paper, the duration of the ki- 
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nation period depends crucially on the parameters of the 
hybrid-like potential, g, m and A-^, as well as on the dis- 
sipation parameter ii. Namely, if g and m decrease (A^ 
increases), reheating becomes less efficient, with the con- 
sequence that the kination period becomes longer. In 
this case, the peak in the high-frequency region of the 
gravitational-wave spectrum becomes more pronounced 
(see Figs, ini [5] and E]). Similarly, if /i increases, the du- 
ration of the kination period also increases, implying a 
higher peak in the spectrum (see Fig. [TU)) . 

In the original quintessential inflationary model 0, 
the kination period is quite long, implying a pronounced 
peak in the high-frequency region of the gravitational- 
wave spectrum [llj. In order to avoid a conflict with 
the integral bound imposed by primordial nucleosynthe- 
sis [sec Eq. (HH)], in the original model one has to require 
that the number of scalar fields reheating the Universe 
is greater than about 10^ meaning that the mini- 
mal Grand Unification Theory is not enough to accom- 
modate quintessential inflation. This situation contrasts 
with the one occurring in the hybrid quintessential in- 
fiationary model. As we have shown in this paper, for 
the allowed values of the parameters of the model, the 
primordial-nucleosynthesis bound is satisfied, by far, by 
any gravitational-wave spectrum. Indeed, even when all 
parameters of the model are pushed to their extreme val- 
ues, the gravitational- wave spectral energy density pa- 
rameter ilcw does not grow above 10~^^. This impor- 
tant difference between the gravitational-wave spectrum 
of the original model and its hybrid variant is due to 
the fact that in the latter the duration of the kination 



period is much shorter, implying a smaller peak in the 
gravitational- wave spectrum. 

The values of the different parameters also determine 
the frequency at which the gravitational-wave spectral 
energy density parameter flew starts to rise at the high- 
frequency region of the spectrum. For extreme values of 
the parameters, this frequency can be as low as 10^ Hz. 
Despite the fact that this frequency falls within the fre- 
quency range at which terrestrial laser-interferometer 
gravitational-wave detectors operate, a detection in a 
near future is improbable, due to the very low value of 
flew, about 10-15 - 10-^^. 

In the gravitational-wave energy spectrum, the maxi- 
mum of the peak is firmly located in the MHz-GHz region 
and, within the hybrid quintessential infiationary model, 
can reach values of about Oqvv — W^^^- A search for the 
primordial gravitational-wave background at this range 
of frequencies was initiated a few years ago, with both 
microwave-cavity detectors jssl l39l | and interferometric 
detectors [i^l . We deem the development of such detec- 
tors to be of the utmost importance, since they will open 
a direct window to the very early Universe, allowing us 
to test our theoretical ideas about the inflationary and 
post-inflationary epochs. 
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